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A microscopic mean-field theory of tfie pliase coexistence between ferromagnetism and supercon- 
ductivity in the weakly ferromagnetic itinerant electron system is constructed, while incorporating 
a realistic mechanism for superconducting pairing due to the exchange of critical spin fluctuations. 
The self-consistent solution of the resulting equations determines the superconducting transition 
temperature which is shown to depend strongly on the exchange splitting. The effect of phase 
crossover from isotropic (Heisenberg-like) to uniaxial (Ising-like) spin fluctuations near the quan- 
tum phase transition is analysed and the generic phase diagram is obtained. This scenario is then 
applied to the case of itinerant ferromagnet ZrZn2, which sheds light on the proposed phase diagram 
of this compound. Possible explanation of superconductivity in UGe2 is also discussed. 
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There has been extensive experimental research done 
recently on the possible coexistence of superconductiv- 
ity (SC) and ferromagnetism (FM) in strongly correlated 
electron materials. First, superconductivity was discov- 
ered in the ferromagnetic metal UGe2 at high pressure 0. 
Later a low-temperature SC phase was found in another 
f-electron compound URhGe[2| and in the d-electron fer- 
romagnet ZrZn2 01 ■ In the best studied cases of UGe2 
and URhGe the experiment strongly suggests that su- 
perconductivity coexists with itinerant ferromagnetism 
in these compounds. Notably, it is the same electrons 
that are involved in both SC and FM, which leads to 
inter-dependence of the corresponding order parameters. 

The aim of this paper is to construct a mean-field the- 
ory of the phase coexistence between FM and SC on the 
border of magnetism, while adopting a realistic mecha- 
nism for superconducting pairing. We show that inter- 
play between FM and SC order parameters has crucial 
effect on the resulting phase diagram. We then propose 
a mechanism that would explain the enhancement of SC 
transition temperature in the FM phase and discuss the 
application of this mechanism to ZrZn2 and UGe2. 

The fact that SC is observed inside the FM region im- 
poses strict limitations on the nature of the SC state. The 
very large internal molecular field due to the exchange in- 
teraction (measured Q to be ^ 240 T in UGe2) excludes, 
due to the Pauli limitation, not only any singlet-pairing 
SC but also any unitary triplet states . In this paper we 
analyse consequences of the so-called non-unitary triplet 
SC state ^ 1^] on the resulting phase diagram, and then 
develop a microscopic theory based on spin-fluctuation 
mediated pairing to proceed beyond thephenomenologi- 
cal level of treatment reported in Ref. |5(. 

A nonunitary triplet state is described Q by the or- 
der parameter AQ/3(k) = (ck,aCk,/3) = [«(d(k) • (T)ay\ap, 
where cr = xcx + ycTj, + iaz denote the usual Pauli 
matrices and the basis of symmetric matrices icray was 



used to represent odd angular momentum pairing. The 
three-dimensional complex vector (i(k) fully character- 
izes the triplet pairing state 0. In what follows we as- 
sume for simplicity an easy axis of magnetization in the 
z-direction. Because of the pair-breaking effect of strong 
exchange field M, only the Cooper pairs with parallel 
spins will survive. In this case of equal-spin pairing we 
can write vector d in the form d = [dx, dy,0). Denoting 
A± = dx i idy, the SC order parameter becomes 
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We shall start from the effective Heisenberg model for 
itinerant electrons with spin s(r) = J^ap V'l(r)cQ^'0/9(r), 
where V'l(r), are electron field operators. Some 

attractive pair-forming interaction V is also assumed: 

HpM+sc = X] 4,QCk,Q - I J dr s(r) • s(r) 

+ \ ^"/3.v(k, k') clka 4/3 Ck'A C-kV- (2) 
k,k' 

Making use of the Hubbard-Stratonovich transformation 
and integrating out fermionic degrees of freedom in order 
to arrive at the effective action in terms of the bosonic 
field operators A±(k) and M(r) = Sz(r), we then deduce 
the mean-field equations for the order parameters in the 
saddle-point approximation. The resulting equations for 
the SC order parameter have usual BCS-like form: 



A+(k) = -i^V^(k,k')i^|g§^A+(k') 



(3) 



where f{E) is the Fermi-Dirac distribution function. 

The magnetic order parameter M enters above 
equations via the quasiparticle spectrum £'±(k) = 



2 



^(ei^ ^ MY + |A±(k)|2. The equation for M looks as 



2M 
where 
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ek ± M . In the hmit of pure magnetism 
(i.e. A± ^ when T 0), this equation reduces 
to the Stoner criterion for itinerant ferromagnet 1// = 

i/^Ek(-|^)-^(o). 

In order to illustrate the interplay between the FM 
and SC order parameters that follows from this model, 
we solved Eqs. (|3I4|I self-consistcntly for the simple case 
of spherical Fermi surface at half filling, while assum- 
ing that SC pairing strength in the p-channel V (k, k') = 



^=i(fc,fc') 



1 

E 



ri™(k)Yi*„(k') has BCS-like form (i.e. 



Vi{k,k') vanishes everywhere except the narrow region 
near the Fermi surface) and does not depend on exchange 
interaction /. The resulting SC transition temperature 
was calculated in the weak-coupling BCS approximation 
and is shown in Fig. ^ as a function of dimcnsionless 
interaction constant / = N(0)I. It is apparent that ex- 
change splitting has large effect on superconductivity, en- 
hancing Tsc in the FM phase for the majority spin chan- 
nel and suppressing it for the minority spin. This is not 
surprising since exchange splitting enhances (suppresses) 
the density of states (DOS) in the majority (minority) 
spin channel, which enters the expression for the dimcn- 
sionless pairing strength Aj = Afcr(0)|Vi|. We note that 
at this stage the symmetry of the magnetic state has no 
effect on Tsc since the mass renormalization effects have 
not been taken into account when calculating Tsc- These 
effects will prove to be very important in what follows. 
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FIG. 1: The calculated SC transition temperature as a func- 
tion of interaction strength / obtained as a result of solving 
self-consistent Eqs. I|8I4^ in the weak-coupling BCS approx- 
imation Tsc ^ 1.14: ujc exp( — 1/Af), where frequency cutoff 
UJc was chosen arbitrarily loc = 0.01 Tf. The solid (dashed) 
line shows Tsc for majority (minority) spin in units of Ef- 
Exchange splitting M is plotted with dotted line (right scale). 

We shall now address the issue of the mechanism of the 
superconducting pairing that arises in the vicinity of the 
magnetic quantum phase transition. It has been shown 
that exchange of spin fluctuations, called paramagnons. 



can lead to an attractive pairing interaction. On the 
paramagnetic side the strength of this interaction was 
derived in the context of superfluidity in '^He 0,11, 13- On 
the FM side of the transition the corresponding formulae 
were obtained in Ref . [lol ITU . Following their approach, 
the (attractive) equal-spin pairing interaction is given by 



'Xa 



1 - I'x?>Xo' 



(-^^T,!), (5) 



where Xo il) = Xo"^ i'i'^i i q) is the Lindhard function of the 
non-interacting system in the given spin channel. The 
usual BCS pairing parameter A in the triplet channel is 
given by A^^^ = N^{0)\V{"^\, where N„{0) is the DOS at 
Fermi level and Vi°'°' is the strength of the interaction in 
the /-orbital channel, which for spherical Fermi surface is 



2kf 



Pi l- 



9k^ 



xV^iuj = 0,q) 
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|k| = |k+q|=fcFa 



where Pi{z) denotes the Legendre polynomial of order 
I. The dominant contribution to the integral comes from 
the small-g region since V'^'^{q) is strongly peaked for q 
0. This allows us to employ the small-g approximation for 
the Lindhard function Xo(q) ^ ^a(0)[l - j^(g/fcF^)^ + 
0(g'')], which enters Eq. |5l 

The mass enhancement near magnetic instability 
renormalizes the BCS pairing parameter Ai to the new 
value A'"* = Xi/Z<t{0) = A^/(1 -h AfJ), where Aq is the s- 
wave pairing interaction parameter and Z(0) is the mass 
enhancement factor j2l| at the Fermi surface given by 



^ ^ ZiO) 
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dY.{Lu, |k| ^ kp] 
du! 



(7) 



where S is the single-particle self-energy. 

It is remarkable that experimentally no SC is observed 
in the paramagnetic (PM) phase of UGe2 and ZrZn2. 
Admittedly, this is different from the result of Fay and 
Appel 01 who obtained comparable values of Tgc on 
both sides of the magnetic transition. Here we propose 
a tentative generic explanation for the strong suppres- 
sion of SC on the PM side of the transition in the nearly 
ferromagnetic metal, which is due to the different na- 
ture of spin fluctuations on the two sides of the mag- 
netic quantum phase transition. We consider a nearly 
ferromagnetic metal which has no preferred magnetiza- 
tion axis in the PM phase and thus is characterized by 
spin fluctuations that are of Heisenberg type. In this 
case both longitudinal and transverse spin fluctuations 
contribute to the effective mass enhancement, so that 
Z{0) = 1 + Xq +Xq . Consequently, the renormalized pair- 
ing strength A* = Xi/Z{Q) and the SC transition temper- 
ature are both small on the PM side. Remarkably, higher 
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FIG. 2: a), b): Calculated properties as functions of 
the exchange splitting I of the itinerant FM with proposed 
Heisenberg- to Ising-type crossover occurring at 7 = 1. Spin 
fluctuations of Heisenberg type are assumed in the PM phase 
and of Ising-type in the FM phase (solid line). This is to 
be compared with purely Heisenberg-like behaviour on both 
sides (dotted line), a) Mass enhancement parameter Za{0), 
where formulae of Ref. llol| were adopted for calculating the 
mass enhancement contributions Aq' and Aq' . b) The effec- 
tive pairing parameter A"^* = X^'^ /Za-{0). c) The calculated 
superconducting Tsc (solid curve, left scale). The frequency 
cutoff LJc, shown in the inset, was extracted from the strong- 
coupling calculation by Levin and Vails [13 and is symmetric 
about 1 = 1 point (Tp is the Fermi temperature). The broken 
curve (left scale) shows Tsc if one adopts the phenomenolog- 
ical formula for ujc from Ref. ITU , and the dotted curve (right 
scale) shows Tsc in the PM phase scaled by a factor of 25. 



values of Tsc on the FM side could be achieved if spin 
fluctuations there were of Ising-type, so that only longi- 
tudinal spin fluctuations contribute to the effective mass 
enhancement Z(0) = 1 + Ag . This would lead to Z{0) 
being about three times larger in the PM phase than it is 
in the FM phase, as illustrated in Fig.[2Ji). The resulting 
pairing parameter A'^* is shown in Fig.jSJs)- We will first 
outline the consequences of this hypothesis on Tsc, and 
then turn to the discussion of its validity. 

The SC transition temperature is notoriously difficult 
to calculate. For the purpose of comparison of Tsc on 
both sides of the FM transition a simple McMillan-type 
formula should suffice: 



Tgc ~ 1.14lJc cxp[-l/A!^' 



(8) 



where cutoff lOc simulates in a crude way the fact that in 
reality l^'^'^(u;,q) is strongly frequency-dependent, being 
sharply peaked at small energy transfers. It turns out 
|l2j that ujc depends strongly on exchange interaction, as 
shown in the inset of Fig. The resulting dependence 



Tsc{I) is plotted in Fig. 13;), which indicates clearly that 
the SC transition temperature is an order of magnitude 
higher in the FM phase than it is in the PM phase. 

Our calculations suggest that Tgc goes through a max- 
imum and then approaches zero at the quantum critical 
point, in accordance with Ref. [llLIT^ . A recent strong- 
coupling calculation by Roussev and Millis 0] suggests 
however that Tsc > generically at the magnetic critical 
point, contrary to our result. We note that though inter- 
esting from the fundamental point of view, the behaviour 
of Tsc directly at the magnetic phase transition is not 
so important in practice, since experimentally magnetic 
transition proves to be first order P, Q , thereby elimi- 
nating the low values of (/— 1) from consideration. 

For outlined scenario to take place, two crucial condi- 
tions are necessary. Firstly, the contribution of soft trans- 
verse spin fluctuations to Aq must be quenched on the FM 
side. This is achieved due to spin waves taking over the 
available phase space as the magnetization M increases. 
Indeed, the fraction of the momentum space available to 
gapless spin fluctuations is q > q± = kpi — kp] , the 
rest being taken by spin waves < q < q±. Long- 
wavelength spin waves themselves do not contribute to 
the mass enhancement in the leading order of the per- 
turbation theory [isj . Thus as the exchange splitting in- 
creases, the soft spin fluctuations shift to larger g- values, 
thereby decreasing their contribution to Aq. 

However, this suppression also affects the pairing 
strength Ai cx x\\ due to longitudinal spin fluctuations 
which become quenched as well. The situation can be 
cured by the second condition: the existence of the quan- 
tum meta-magnetic transition (MMT) jl6| somewhere in 
the FM phase. Indeed, the longitudinal susceptibility 
X|| is peaked near the jump in magnetization accompa- 
nying such a transition, as seen experimentally in UGe2 
|17| and in Sr3Ru2 07 [l3- As a result, longitudinal spin 
fluctuations will be enhanced and the material will ap- 
pear effectively Ising-like near the MMT, justifying the 
assumption made above. However unlikely the "coinci- 
dental" presence of meta-magnetic transition near the 
quantum transition to the FM phase may appear at first 
sight, the experiment suggests that this is not uncommon 
in the ferromagnetic strongly correlated electron materi- 
als. Indeed, the MMT has been observed in Sr3Ru207 
El, UGe2 [H 113, and recently in ZrZnz 

We now turn to the application of the above model 
to the experimentally studied materials. Figure |21 shows 
the generic phase diagram of an itinerant ferromagnet 
that arises from studies of ZrZn2 0, 0, Ull and UGe2 
P, ll^, HE 13 • The Curie temperature Tc is suppressed 
to zero at pressure pc, where the transition appears to be 
first order [l],|3{ in both compounds. Another feature, the 
crossover line between the two ferromagnetic phases, 
FMl and FM2, is also shown. This crossover exhibits 
itself as an anomaly in the measurements of resistivity 
Hm m and specific heat ,26, 211 in UGe2, where it 
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FIG. 3: The schematic p-T phase diagram of a generic itiner- 
ant ferromagnet that exhibits a meta-magnetic transition at 
p = Px, associated with the SC region inside the FM phase, 
based on studies of ZrZn2 and UGe2. Pressure px is negative 
in ZrZn2 and positive in UGe2 . The dependence of low tem- 
perature magnetization, M, on pressure is shown schemati- 
cally with dotted line I2a.l2al . Tc is the Curie temperature; Tx 
is a crossover line between two different ferromagnetic phases, 
FMl and FM2. The filled circles denote critical endpoints, 
below which the transitions are first order. The dashed line 
indicates a crossover, rather than a sharp phase transition. 



occurs at pressure Px ~ 12 kbar. In the case of ZrZn2 Px 
appears to be negative 0,112 and thus cannot be probed 
directly, however the MMT can be rendered to positive 
pressures by applying external magnetic field, where it 
has been studied in the de Haas-van Alphen [Ul and 
magnetization measurements |22j . 

It is notable that SC is observed only on the FM side of 
the transition in ZrZn2 Q , but not in the PM phase. This 
fact can be readily explained by the scenario proposed in 
this work. Indeed, ZrZn2 is a three-dimensional ferro- 
magnet with cubic symmetry, from which the isotropic 
(Heisenberg) nature of spin fluctuations in the PM phase 
is deduced. By contrast, it follows from the proposed 
model that spin fluctuations develop Ising-like symm etry 
upon entering the FM phase, where the existence |2lll22l| 
of meta-magnetic transition ai p ~ Px is crucial to our 
argument. It thus follows from our calculations that the 
SC transition temperature must be strongly suppressed 
in the PM phase of ZrZn2. 

The situation is more intricate in UGe2, where evi- 
dence of strong uniaxial anisotropy exists on both sides 
of magnetic transition at Pc [ll 1131 • However in the light 
of recent measurements of specific heat it becomes 
evident that the very narrow SC region is centred around 
the MMT at px rather than pc- This can be easily un- 
derstood given that the transition at pc is strongly first 
order and therefore has a strong pair-breaking effect 
on SC. It is hence not surprising that no SC is seen on 
both sides of Pc- The presence of SC in UGe2 is instead 
due to critical spin fluctuations at px, which is only a 
weakly flrst order transition. Wc note that this view has 
been expressed already in earlier works on UGe2 [23ll2^ . 
Our proposed theoretical model has thus an indirect ap- 



plication to UGe2 in a sense that the observed uniaxial 
(Ising) symmetry enhances TgCi which would have been 
much more suppressed if the Heisenberg-like spin fluctu- 
ations had prevailed in this compound. 

The above argument already suggests that SC phase 
must be suppressed both in FMl and PM phases of UGe2 
close to Pc- It should be noted that the qualitative change 
in the Fermi surface observed at the magnetic transition 
by dc Haas-van Alphen experiment may be another 
factor that suppresses SC near pc- In this context, the 
existence of a double peak structure in the electronic den- 
sity of states very close to the Fermi level has been pro- 
posed ll^l as a possible microscopic explanation. We also 
note that absence of SC in the PM phase of UGe2 may 
be partly due to the spin degeneracy of the Fermi sur- 
face as it can, in principle, enhance spin-flip processes of 
the electrons forming a Cooper pair, which would have 
detrimental consequence on spin-triplet SC state. 

In conclusion, we have formulated a mean-field the- 
ory of coexisting FM and SC in terms of the equations 
for the corresponding order parameters that have to be 
solved sclf-consistently. We have also incorporated a mi- 
croscopic mechanism of the SC pairing due to the ex- 
change of spin fluctuations in our model. A scenario 
based on Heisenberg- to Ising-type crossover has been 
proposed, which provides a natural explanation of the en- 
hancement of SC on the FM side of magnetic transition, 
observed experimentally in ZrZn2. The apparent sup- 
pression of SC in the PM phase of UGe2 is explained by 
the detrimental effect that the strongly first-order phase 
transition at pc has on pair-forming spin fiuctuations. 
The proposed theoretical model supports the evidence 
of SC in UGe2, as superconductivity is predicted to be 
enhanced by Ising-like spin fiuctuations near px- 
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